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Abstract 

The standard block orthogonal (SBO) polynomials Pi-n{x), < i < n are real polynomials of 
degree n which arc orthogonal with respect to a first Euclidean scalar product to polynomials of 
degree less than i. In addition, they are mutually orthogonal with respect to a second Euclidean 
scalar product. Applying the general results obtained in a previous paper, we determine and 
investigate these polynomials when the first scalar product corresponds to Hermite (resp. Laguerre) 
polynomials. These new sets of polynomials, we call Hermite (resp. Laguerre) SBO polynomials, 
provide a basis of functional spaces well-suited for some applications requiring to take into account 
special linear constraints which can be recast into an Euclidean orthogonality relation. 

PACS numbers: 02.10.Ud, 02.30.Gp, 02.30.Mv, 21.60.-n, 31.15.Ew, 71.15.Mb 
1. Introduction 

In a first paper [1] (paper I), within the general frame of an Euclidean vector space endowed with distinct 
scalar products, and using the Gram-Schmidt orthogonalization (G-SO) process, we investigate the general 
properties of standard block orthogonal (SBO) polynomials. By definition, for a given non-negative integer i, 
these real polynomials Pi-n{x), n ^ i,i + I, . . of degree n, are 

(i) orthogonal to the i-dimensional subspace Vi of real polynomials of degree less than i with respect to a 
first Euclidean scalar product ( , ) defined by a non-negative weight function w{x) on the interval [a, 6], 

(x"* , Pi;„) := / dx w{x)x"' P^,n{x) = m = 0,...,i-l n = i,i + 1, . . . ; (1.1) 



(ii) mutually orthogonal with respect to a second Euclidean scalar product (, )2, similarly defined by a 
distinct non-negative weight function W2 (x) on the same interval ^ , 

{Pi;ni , Pi;n)2 := / dx W2{x)Pi,jn{x)Pi-n{x) = Hi-nSm^n m, n = i , i + 1 , . . . . (1.2) 
J a 

In other words, for a given subspace Vi, there exists a unique orthogonal complement of codimension 
i, orthogonal with respect to the scalar product ( , ), and complementary with respect to the space of real 
polynomials V of any degree ^, 

{V^,Vi^)=0 and V^V^®Vt- (1-3) 

Then, {Pi-^i, Pi-i+i, ■ ■ ■} is an orthogonal basis of V^^ with respect to the scalar product ( , )2. These basis 
polynomials are defined uniquely apart from a constant multiplicative factor for each polynomial. 

^ More generally, the scalar product ( , ) can be defined by the positive measure {D, d/i} where, is a 
non-constant and non-decreasing real function on the real domain V. Similarly, ( , )2 can be defined by 
{I?2, c?M2}- See section 3.2 in paper I. 

^ See equation (3.4) in paper I, written with Vn instead of V and valid for any non-negative integer value 
of A^. 



Classical orthogonal polynomials (i.e., the Hcnnitc, Lagucrrc and Jacobi polynomials) are orthogonal 
with respect to a scalar product defined with a particular positive measure. These polynomials arise fre- 
quently and have been studied in great detail [2-5] ^. The purpose of this paper is to determine and 
investigate the SBO polynomials in the special cases where the first measure, {[a, 6], w{x)dx} corresponds to 
the Hermite and then the Laguerre polynomials. This work has been motivated by the recent consideration 
of ^constrained orthogonal polynomials' by Giraud et al [6-8]. Their physical motivations, sketched in section 
1 of paper I, were basically to define orthogonal basis fimctions, taking into account special linear constraints 
(i.e. orthogonality relations), to provide a functional space well-suited for some specific applications. Thus, 
aiming to use the Hohenberg-Kohn variatonal principal [9] for the ground state energy in nuclear physics, 
the variable function, i.e. the particle density, has to satisfied the constraint of particle number conservation. 
In toy models, expanding the particle density on SBO polynomials, these authors only consider the case 
i = 1 (and partly i = 2) for the Hermite polynomials in [6, 8] and the Laguerre polynomials in [7]. Our 
works generalize these results to all non-negative integer values of i and provide a systematic study of the 
properties of these polynomials, in particular getting new recurrence and differentiation formulae. 

This paper is organized as follows. Explicit expressions and properties of the Hermite (resp. Laguerre) 
SBO polynomials are given in section 2 (resp. 3), i.e. for: the connection coefficients in section 2.2 (resp. 
3.2), the difi'erentiation formulae in section 2.3 (resp. 3.3), the recurrence relations with respect to the degree 
n in section 2.4 (resp. 3.4) and with respect to both i and n in section 2.5 (resp. 3.5), the special value 
at a; = in section 2.6 (resp. 3.6), the properties of the zeros in section 2.7 (resp. 3.7), and finally, the 
interrelation between the Hermite and the Laguerre cases in section 3.8. We give our conclusions in section 
4. For completeness and to illustrate how the general methods used allow to recover well-known results in a 
unified way, several properties of the Hermite (resp. Laguerre) polynomials are briefly recalled in appendix A 
(resp. B). This aims also at emphasizing the similarities to and the diS'erences from the study of the classical 
orthogonal polynomials and the SBO polynomials. Some expressions related to the Pochammer symbol and 
generalized hypergeometric series arc computed in appendix C. Some determinants with entries expressed 
in terms of gamma function, except possibly for the last row, are evaluated in appendix D. Finally, some 
explicit special cases for the Hermite (resp. Laguerre) SBO polynomials are listed in appendix E (resp. F). 

Throughout the remaining of this paper, the same conventions and notations as in paper I are used 
in addition 

- references to (sub) sections, appendices and equations of paper I are preceded by I-; 

- [x] denotes the largest integer less than or equal to the real number x; 

- for any function f : x ^ f{x), the first (resp. second) derivative with respect to x is denoted by a prime, 
/' (resp. a double prime, /"). 

2. Hermite standard block orthogonal polynomials 

Definition and main properties of Hermite polynomials, Hn, n = 0, 1, . . ., are recalled in appendix A. They 
correspond to the measure defined by 

o := -oo 6 := DO w := e""^^ . (2.1) 



2.1. Metric tensor components for a particular weight function W2 
(i) Basic quantity ^j,k- With the same interval, [a, 6], let us choose, 

W2 = wi^^ := e""^' n>0. (2.2) 

^ The standard textbook is [2]. See, also [3] chapter 22, [4] sections 10.6-10.13 or [5] section 8.9. 

^ See, [1] at the end of section 1. In particular, it is recalled that: - monic polynomials, i.e. with the 
coefficient of jhe highest power equal to one, and also any related quantities are denoted by hatted letters, 
e.g., P„ and hn, - classical orthogonal polynomials are defined according to [2] chapter 22, [3] chapter X or 
[4] section 8.9. 
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This choice is motivated by the apphcation considered in [6-8] for the special value //, = 2. Our study will 
be restricted in (ii) below to this special case, since it allows to push the calciilation to its end analytically. 

As shown in section 1-3.3.3, w and W2 being even with a = —b, the SBO polynomial Pi;„, n = + . . ., 
as the Hermite polynomials, are even or odd according as n is even or odd, equation I-(3.41), 

PrA-x) = {-irPi-n{x). (2.3) 

Actually, the even and odd polynomials separate from each other and one has, equation I-(3.50), 

if (-1)*+" = 1 then Pi_i;„ = Pi;„ < i - 1 < n . (2.4) 
From equation I-(3.12), the basic quantity Jj,k{l^) reads 

/CXD 
dxe-''='"Hj{x)Hk{x) (2.5) 
-oo 

which vanishes by parity if j + k is odd. Expanding the Hermite polynomials in terms of powers of x with 
the coefficients given by equation (A. 7) and using equation (A. 4), one finds 

72i,2fe(M) = (-i)^+'=22W+'=)r(i + i)r(i + k) 



X 



3 /A ^^ly+g 



^oi^W Wr(i+p)r(i + g) i_ 



= (_i)i+'=220+fc)r(| + i)r(i + 



j k 

X 



The value of the sum over p and q above is considered in lemma C.2 (appendix C). It can be expressed in 
terms of an hypergeometric function or the function S defined by equation (C.8), then, 

l2jM^^) = 22(^+'=) r(i + j) r(i + fc) m"* s{j, k, ^■-^^-\-^,-^) . (2.7) 

For the odd polynomials, from equation (A. 8), a similar calculation yields 

72,+i,2fe+i(M) = (-1)^+^- 22(wi) r{i + j) r(f + k) fi-i s{j, fc, I; . (2.8) 

Alternately, these last two equations can be obtained using, with some care, the definite integrals given in 
[5, 10] ^. Another way to get the basic quantity 7j,fc(/x) is to compute the scalar product (F, , Ft)2 where, 
Fz{x) is the generating function of Hermite polynomials given by equation (A. 29). Then, 'yj,k{lj) is related 
to the coefficient of sH'' in the expansion of this scalar product, i.e., 

/oo ^1 
dx e-'^-" F,{x)Ft{x) = J2 -^^^JA^^>'t'' 

= n%xp?£i±il^i^li£^. (2.9) 

^ See, [5] equation 7.374 5 where, the second argument of the hypergeometric function should be — n 
instead of n (the result has to be symmetric in m and n) and in addition for real, a has to be taken 
positive. See also [10] 2.20.16 1 where, the factor (±1)™+" should be replaced by (-1)™+", or else, 2.20.16 
6. In both cases, one has to use also the linear transformation formula [3] 15.3.6. 



For /i = 1, i.e. for u>2 = w, one has to recover the orthogonahty condition of Hermite polynomials. 
Indeed, from equation (C.14), 



S{j, k, i; -1, -1) = YfJ^S,,k (2.10) 
whence, with the duplication formula [3] ^, 

72,,2fc(l) = 24Ji!r(i + j) S,,k = 2^^{2jy.nhj,k (2.11) 

recovering thereby the value of h2j given by equation (A. 6). The calculation is similar for 72^+1, 2fc+i(l)- The 
results follows readily also from equation (2.9), indeed, {Fg , ^4)2 = Tr^e^** = YlJLo ^ji^'^Y / U^-)^ ■ 

(a) The special case W2 = ^2^'' = e~^^^. When ji = 2, one gets simple expressions. Indeed, from equation 
(C.15) setting z = — one has 

72i,2fe(2) = (-1)^+'= 2^+'=-^ r(i + i + fc) (2.12) 
72,+i,2fc+i(2) = (-l)J-+'=2J-+*=+4r(|+i + fc) (2.13) 



or, equivalently [5] ^, 



^.^^{2) = { (-l)[^/^™2iO+'=-i)r(i(j+fc + l)) ifj + fceven ^^.M) 
otherwise. 



Alternately, this result follows also from equation (2.9), 

00 2n 



« . ^.)^ - G) = (I) ' i: i: (-'r-^- „i„ffL)i ^ ("^) 

Setting j = m, k = 2n — m and using the duplication formula to get (2n)!/n! = 7r~5 2^"r(i + n), one recovers 
equation (2.14). 

From now on, except explicitly noted, only this special case pL — 2 will he considered. 

Then, the corresponding polynomials Pi-^m < i < n are called the Hermite SBO polynomials. The case 
i = 1 has already been considered in [6, 8] with different choices of weight functions and normalization. (See 
footnote 50 in appendix E for interrelations with our polynomial -Pi;™-) The case i = 2, only for the odd 
polynomials, has been also considered in [6], as already discussed in the conclusion of paper I (see footnote 
44 of paper I). (See footnote 51 for interrelations with our polynomial P2;2n+i.) 

2.2. Connection coefficients between block orthogonal and classical polynomials 

The relevant quantities are given by equations I-(3.42)-I-(3.47) in terms of determinants with the {j, k)- 
element of the type r{c + j + k), except possibly for the last row. This kind of determinant is evaluated in 
lemma D (appendix D). From equations I-(3.43), I-(3.46), (2.12), (2.13) and (D.3), one finds for < i < n, 

n 

^e) ^ _ r(i + z + i) Zg_,:=l (2.16) 

71 

ZS = 2(^+"+^)("-+i) Hij-iy.ril + i+j) Z^_, := 1 (2.17) 
i=i 

6 See, e.g., [3] 6.1.18, T{z)T{\ + z) = 2-2^+i7r5r(2^) . 

^ Equation (2.12) is given in [5] equation 7.374 2, with the wrong sign: the result should be symmetric in 
j and k and positive for j = k. 



and therefore, from equations I-(3.44) and I-(3.47) with A;„ given by equation (A. 3) [6] ^, 

^2i-l;2n = ^2i;2„ = 2-(2"+i)(n - i)\ T{^+i + n) (2.18) 

^2i;2n+i =i?2m;2«+i =2-(2"+i)(n-i)!r(|+i + n). (2.19) 
The monic SBO polynomials follow from equations I-(3.42), I-(3.45) and (D.2), for < i < n, 

P2.-l;2„ = P2,2„ = E 2-("+™) f " - ') Ilf±l±4ff2™ (2.20) 
P2i;2n+1 = P2i+l;2n+l = E 2-(l+"+™) {^^ ^^'^ ^ (2.21) 

where, according to equations I-(3.13), (2.4) and (A. 3), for i > 0, 

Pi,i = Hi := 2-' Hi Pi,i+i = Hi+i = 2-('+i) Hi+i . (2.22) 

The orthogonality and normalization of these polynomials with respect to the scalar product ( , )2 can be 
checked directly using equation (C.14), e.g., for the even polynomials with equations (2.20) and (2.14), 



/. 



oo 

2x 



dxe P2i-2n{x)P2i:2n'{x) 

OO 

= 2-("+"'+5)r(i +i + n)T{\+i + n')S{n -i,n' 2i; -1; -1) 

= 2-(2"+5)(n - i)\ r(i + i + n)6n,n' (2.23) 

recovering equation (2.18). 

Equations (2.20) and (2.21) can be inverted to give for < i < n, 

H2n = E(-l)"-'"2"+-f"-^.)lfi±^P.,2„. (2.24) 

H..» = t(-l).-™2'«-(7_-)i|±l±^P„«,„«. (2.25) 

As for the Hermite polynomials in appendix A, these relations can be obtained either using the general 
formulae I-(A.22) and I-(A.28) where, the determinants which occur can be evaluated with equation (D.2), 
or directly, using lemma A (appendix A), e.g., for the even polynomials, setting j = m — i, k = n — i, 
pj = 2-(*+J)r(i + 2i+j) and aj = (-1)^'2*+Jr(i + 2i + j). 

It should be noted that the connection coefficients relating the polynomials Pi-n and Hn in equations 
(2.20), (2.21), (2.24) and (2.25) have the same structure as the connection coefficients relating iJ„ and the 
monomials in equations (A. 7), (A. 8), (A. 12) and (A. 13), respectively. This is due to the fact that in both 
cases, the metric tensor components are of the same type T{c + j + k), apart from some factors depending 
separately on the indices j and k, see equations (2.14) and (A. 4). (See also equation (2.30) below.) 

As already noted in section 3.2 of paper I (see equations I-(3.18) and I-(3.19)), for i = 0, Po-n has to 
be the unique standard monic orthogonal polynomial Q2-n corresponding to the weight function W2. Since 

'W^\x) = 'w{\/2x), with the normalization (A. 3) and changing x into \/2x in the integral defining the scalar 
product of Hermite polynomials, one gets for n > 0, 

Po,n{x) = 2-5"#„(V2x) = 2-i'^Hn{V2x) . (2.26) 

* These equations generalize the formulae given in [6] (4) (resp. (19)) for the special case i = 1 (resp. 
i = 2 for the odd polynomials) with different choices of w, W2 and normalization, see equation (E.12) and 
footnote 50 (resp. 51) 
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Thereby, sotting « = in equations (2.20), (2.21), (2.24) and (2.25), using the duplication formula and 
changing m into n — m in the sum over m yield the connection coefficients between the two basic sets of 
polynomials {Hm, m = 0, . . . , n} and {HmiV^x), m = 0, . . . , n} for n > 0, 

[n/2] , 

H4V2x) = 2^("-^"-) ^. Hr,_2m{x) (2.27) 

^ ' ^ ml in — 2m)! 

m=0 ^ ' 

["/2] , 

iin{x) = (-1)'"2-^" g„-2„^(^/2a;) . (2.28) 

m!(n — 2m)! ^ ' 

m=0 ^ ' 

These rcmaikablc; relations can be checked using the integral [5, 10] ®and the orthogonality relation of 
Hermite polynomials with /i„ given by equation (A. 6), e.g., one has 

r dxe--'HniV2x)H„_2mix) = 2^^^-^^)—/^ = 2^("-^") , . h^.2m ■ (2.29) 

J_^ ml m\[n — 2m)\ 

Another way to get the relation (2.27) is as follows. From equations (A. 4) and (2.14), one has 

{x^ , x'') = (-l)y/21+[fe/2]2-i0+'=-i)7^. fe(2) = V2(^{-l)^/^h-i^Hj , {-l)i''/^h-i''Hk) ^ (2.30) 

namely, the metric tensor components of the scalar product ( , ) in the basis {a;", n = 0, 1, . . .} and the 
metric tensor components of the scalar product ( , )2 in the basis {(— l)["/^l2~2"iy^, n = 0, 1, . . .} are equal 
apart from a constant factor. Hence, while G-SO with respect to ( , ) of the basis {a;", n = 0, 1, . . .} defines 
the Hermite polynomial iJ„, G-SO with respect to ( , )2 of the basis {(— l)["/^l2~2"i/„^ n = 0, 1, . . .} defines 
the polynomial Pq-u- In other words, the linear mapping defined by x" i-^ (— l)["/^l2~i"ff„, n = 0, 1, . . ., is 

such that Hn i— > oc Po,n, n = 0, 1, Therefore, apart from a multiplicative factor c, Po,n reads as a linear 

combination of (-l)['"/2l2-5™iy„, m = 0, 1, . . . with the same connection coefficients as in the expansion of 
Hn in terms of the monomials i.e. from equation (A. 9), 

^0;- = ^ E u „ ((-l)'("-^"')/^'2-^(»-^'»)g„_2J . (2.31) 
^ m! n — 2m)!v / 

m=0 ^ ' 

With Po:n a monic polynomial, computing the coefficient of fixes c = (— l)l"/^l2^^"n!, recovering thereby 
equation (2.27) with (2.26). 

From equations (2.4) and (2.26), one has 

A;2n+i(a;) = Po;2„+i(a;) = 2-3("+^)if2„+i(\^a;) (2.32) 

(see also the comments about equation I-(3.52) and in the conclusion of paper I, Pi;2ra+i = Q2;2n+i 

). This 

relation can be checked setting i = in equation (2.21) and comparing with the expansion of i?2n+i(\/2a;) 
given by equation (2.26). 

It follows from equations I-(3.10), I-(3.11), (2.20), (2.21) and (A.IO) that, 

Ri,n = (2.33) 
S2i-i;2n = S2i;2n = " J («(2i - 1) + n(2n - 1)) (2.34) 

S2i;2n+i = S2i+i;2n+i = - J (i(2z + 1) + n(2n + 1)) . (2.35) 

^ See [5] equation 7.374 4 or [10] 2.20.16.7. p 503. 
It is underlined in section 1-2.1 that all the relevant quantities in G-SO procedure depend only on the 
metric tensor components. 
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For 1 = 0, one gets from above 5*0:11 = —ri(ii — l)/8, which ean be checked using equations (2.26) and (A. 10). 

For the general properties of SBO polynomials about the projection operators, the integral representa- 
tions and the normalizations, see sections 3.3.1, 3.3.2 and 3.3.4 in paper I, respectively. 

2.S. Differentiation formulae 

Let us summarize well-known results for all classical orthogonal polynomials; 

• as for any standard orthogonal polynomials, the three term recurrence formula I-(3.55) holds; 

• the weight function which defines the scalar product satisfies a differential equation. This implies several 
differentiation formulae, e.g., the linear homogeneous differential equation of the second order , 

A{x)y" + B{x)y' + Xny = (2.36) 

where, A{x) and B{x) are independent of n and A„ is independent of x. 

All these relations follow from the three term recurrence formula and one of the first order differentiation 
formulae. Moreover, among the three formulae — the recurrence equation, a first order differentiation equation 
and the second order differential equation — any pair of equations implies the third one. In particular, two 
first order differentiation formulae or, the first and second order differentiation formulae, imply the recurrence 
formula. This is the method used to get a five term recurrence formula in the next subsection. However, while 
the differentiation formulae for classical orthogonal polynomials can be derived in two ways, as recalled in 
appendices A (item 4) and B (item 4), the second method only can be used for block orthogonal polynomials. 
This is due to the fact, already encountered in sections 3.3.5 and 3.3.7 of paper I, that the scalar product 
(, )o, defined by equation I- (3. 5) is no longer defined by a single integral. Here, for a given nonzero i, 
one has to deal with both scalar products ( , ) and ( , )2, each of them defined by an integral with a different 
weight function. 

For the scalar product ( , ), the weight function w satisfies the differential equation (A. 17), which implies 
the relations (A. 18) and (A. 19). Similarly, for the scalar product (, )2 associated with the weight function 
^2'*^ defined by equation (2.2), one has for any polynomials / and g, 

wi"^' = -2fixwi''^ (2.37) 
{f,g')2=-if'-2fixf,g)2 (2.38) 
(/" - 2/.X/' , 5)2 = -(/', fl')2 = (/, fl" - 2i,xg')2 . (2.39) 

(i) For < £ < i, the polynomial H'^ of degree — 1 is in the subspace Vi- Therefore, from the orthogonality 
condition (1.1), H'^ and Pi-n for i < n are orthogonal with respect to the scalar product (, ), and from 
equation (A. 18), one has 

= (Pi;„ ,Hi)=- [Pi.^ - 2xPi,n , He) . (2.40) 

Hence, the particular polynomial — 2xPi-n of degree n + 1 is orthogonal to the subspace Vi with respect 

to the scalar product (, ), i.e. it belongs to the orthogonal subspace V:^, see equation (1.3), and thus 
reads ^i;n,jPi;j- Computing the coefficient of a;"+^ in P/.„ — 2xPi-n fixes ai-n,n+i = —2. By parity, 

cti-n,n, c(i-n,n-2, ■ ■ ■ vanish. Now, with respect to the scalar product ( , )2, one has from equations (1.2) and 
(2^38), for'j = i,...,n+l, 

= -{P;.^ - 2(/X - l)xh,, , h,n)2 • (2.41) 

" See, e.g., [4] 10.6 (ii) p 164. 

Any Euclidean scalar product (, )i defined on Vi, together with (, )2 defined on V^- induce a new 
Euclidean structure on V with the scalar product defined by {p , q)o '■= {Pi , qi)i + {pf' , Qi~)2 for any p = 
Pi+pi eV, p^eV^, pi eV^ a.ndq = qi + q^ GV, QiGVi, eV^. 

Actually, equation (2.40) still holds for i = i. Since from equation (2.22), Pi-i = Hi, the component 
0!i;n,i considered in equation (2.41) below, vanishes. This result is obtained also by parity. 
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In the special case /x = 2, it turns out that P/.^- — 2(/i — l)xPi-j is the particular polynomial P/.^- — 2xPi-j, 

introduced above for j = n (i.e. of degree j + 1 in the subspace V^), and which reads J2k^i (^i;j,kPi;k- Thus, 
for II = 2, considered only henceforth, one has for i < j < n + 1, 

^i;n,j-^i;j ~ (^i;j,n-^i;n ^ ^ ^k,n • (2.42) 
k=i 

For n = i, all the a-coefEcients are known: ai-i^i = 0, Q:i;i,i+i = —2. For n > « + 1, it follows from above 
that for n > i + 2 and j = i, . . . ,n — 2, on-n,] vanishes. It remains to consider the case j = n — 1. With 
cti;n-i,n = "2, equations (2.18) and (2.19), one gets (by inspection of the four cases according to the parity 
of i and n), 

Ki;n := a^;n,„-l = 2^^ = " • (2-43) 

Finally, for ;U = 2 and < i < n, one finds the first order differentiation formula 

P/.„ = -2Pi;„+l + 2xPi;„ + Ki;„Pi;„-l (2.44) 

which still holds for n = i, since then Ki-n vanishes. This equation has to be compared with equation (A. 26). 
For n = i, both equations coincide, using equation (2.22). For i = 0, using equation (2.26), one gets 

= —Hn+i + xHn + nHn-i (2.45) 

which follows from, e.g., the differentiation formula (A. 21) and the recurrence formula (A. 16). 

With the boundary value of Pi-i, given in terms of Hermite polynomials by equation (2.22), the differ- 
entiation formula (2.44) can be used to get step by step Pi;„ for all i and n with < i < n. One can also 
take advantage of the parity relation (2.4) to divide by two the number of polynomials to be computed. A 
table of polynomials Pi-n for i = 0, 1, 2 and n = i, i + 1, . . . , i + 4 is given in appendix E. 

Following the same steps as above, for < £ < i < n, the polynomial xfl^ of degree I is orthogonal to 
Pi-n with respect to the scalar product ( , ), and from equation (A. 18), one has 

= {Pi,n , xH[) = [xPi-n ,H'i)=- (xPi.n + Pi;n ' '^x'^Pi.n , He) = - (a;P^.„ - 2a;2p,.„ , Hi) . (2.46) 

Hence, the particular polynomial xP/.„ — 2 x^ Pi-n of degree n + 2 belongs to and reads A;n.j Pi-j- 

Computing the coefficient of in a;P/.„ - 2a;^Pi;„ fixes /?i;„,„+2 = -2. By parity, f3i-,n,n+i, f3i-,n,n-i, ■ ■ ■ 

vanish. Now, with respect to the scalar product ( , )2, one has from equations (1.2) and (2.38), for i < j < 
n+1, 

Pi;n,j Hi- j = (Pi;j , xP^.^ — 2x Pi-n )2 = {xPi-j , Pi-n) 2 ~ {"^X Pi-j , Pi;n) 2 

= - {xPlj - 2(/z - l)x2p,., , P,„)2 - H^,nS,,„ . (2.47) 

In the special case /i = 2, it turns out that xPlj — 2(/i — l)x'^Pi-j is the particular polynomial xPlj — 2x^Pi-j , 

introduced above for j = n (i.e. of degree j + 2 in the subspace P/"), and which reads X);^^^ Pi;j,kPr,k- Thus, 
one has for i < j < n + 1, 

i+2 

f3i-,n,jHt.j = ^— /3i;j,„ ^ (5fe,„ — Sj^n^Hi-^n (2.48) 
k=i 

This generalizes the formula given in [6] (7) for the special case i = 1 with different choices of w and 
W2, see footnote 50. 
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and therefore, = — i. Then, for n = i and n = i + 1, all the /3-coefEcicnts are known: /J^-i^i = 

— = 0, /3i;i,i+2 = —2 and /3i;i+i,i = 0, Pi-i+i^i+i = —5, (ii-i+i,i+2 = 0, /3i;i+i,i+3 = —2. For 

n > i + 2, it follows from equation (2.48) that for n > i + 3 and j = i, . . . ,n — 3, /3i-n,j vanishes. It remains 
to consider the case j = n — 2. Prom Pi-n-2,n = ~2 and equation (2.43), one gets 

Pi;n,n-2 = 2 I^' " = ^Ki-nl^i-n-l ■ (2.49) 
Hi-n-2 ^ 

Finally, in the case ^ = 2, for < i < n and «;,;„ given by equation (2.43), one finds another first order 
differentiation formula, 

XP!.^ = -2Pi;„+2 + ^(4a;2 - l)Pi,„ + \l^i;nKi;n-lPi;n-2 (2.50) 

which still holds for n — i and n — i + 1, since then Ki-nKi;n-i vanishes. This equation, which relates 
polynomials having the same parity, has to be compared with equation (A. 27). For i = 0, using equation 
(2.26), one gets 

AxH'^ = -Hn+2 + 2(2^2 - l)iJ„ + 4n(n - l)i/„_2 (2.51) 

which follows, e.g., from the differentiation formula (A. 27) and the recurrence formula (A. 28). 

Starting with the values of P,;, and Pi-i+i given by equation (2.22), the differentiation formula (2.50) 
allows, as equation (2.44), to compute step by step Pi-n for all i and n with Q <i <n (see appendix E). 
(ii) Following the same arguments again, for < £ < z < n, the polynomial H" — 2xHf of degree i is 
orthogonal to Pi;„ with respect to the scalar product ( , ), and from equation (A. 19), one has 

= (P,.„ , H'l - 2xH't) = (P/.; - 2xi^.„ , He) . (2.52) 

Hence, the particular polynomial P/.'„ — 2a;P/.„ of degree n belongs to P/- and reads J2^=i1i;ri,jPi;j ■ Com- 
puting the coefficient of a;" in P/.'„ — 2a;P/.„ fixes 7i;„,„ = — 2n. By parity, Ji;n,n-i,'li;n,n-3, ■ ■ ■ vanish. With 
respect to the scalar product ( , )2, one has from equations (1.2), (2.38) and (2.39), for i < j <n, 

li;n,j Hi.j = {P,,, , ^ - 2xPlJ^ = (P;, , P;:^ - 2fixPlJ^ + 2{^i - l)(xP;, , P^ 

= ((^^ - 2(m - l)xP^.j) - M^P'i;j - 2(m - l)a;'-Pi;i) . hn)^ - 2(m - mi;n5j,n ■ (2.53) 

In the special case 1-1 = 2, once again, P/'^- — 2(/U — l)xPlj (resp. xPlj — 2(/x — l)x'^Pi-j) is the particular 
polynomial P/'^- — 2xPlj (resp. xP-.j — 2x^Pi-j), introduced above for j = n (i.e. of degree j (resp. j + 2) in 
the subspace P/"), and which reads X^^^j li;j,kPi;k (resp. Y^^i Pi;j,kPi;k)- Thus, one has for i < j <n, 

IvjUj-^i-J ~ {li;j,n ^ ] ^k,n ^Pi;j,n ^ ^ ^k,n '^^j,njHi;n' (2*54) 
k=i k=i 

For n = i and n = z + 1, all the 7-coefficients are known: ji-i^i = —2i and 7i;i+i,i = 0, 7i;i+i,i+i = —2{i + 1). 
For n > 2 + 2, it follows from above that for n > i + 3 and j ^ i, . . . ,n — 3, ^i;n,j vanishes. It remains to 
consider the case j = n — 2. From Pi-n-2,n = —2 and Pi-n,n = one gets 

H- 

7i;n,n-2 = -^0i;n-2,n q = 2Ki;„Ki;„_l . (2.55) 

Hi-n-2 

Finally, in the case /i = 2, for < i < n and Ki-n given by equation (2.43), one finds the second order 
differentiation formula 

- 2a;P^.„ + 2nPi;„ = 2Ki;„Ki;„_i Pi;„_2 (2.56) 

This generalizes the formula given in [6] (15) for the special case i = 1 with different choices of w and 
W2, see footnote 50. 
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which stih holds for n = i and n = i + 1, since then Ki-nKi-n-i vanishes. This equation has to be compared 
with the differential equation (A. 24). For n = i and n = i + 1, both equations coincide, using equation 
(2.22). For i = 0, from equation (2.26), one gets 

H^-xH^ + nHn = 2n{n - 1) Hn-2 (2.57) 

which follows, e.g., from equations (A. 22) and (A. 24). 

2.4- Recurrence relation with respect to the degree n 

In the case /i = 2, eliminating between the two first order differentiation formulae (2.44) and (2.50) 
yields for < i < n the five term linear recurrence formula, 

Pi;n+2 ~ + —Pi-^n + 2 ~ 2 = (2.58) 

where, Ki-n is defined by equation (2.43). This relation still holds for n = i and n = i + since Ki-n vanishes 
for n = i. It has to be compared with the three term recurrence formula (A. 16). For i = 0, with equation 
(2.26), one finds 

if„+2 -2x Hn+i + 2 + 4n a; - 4n(n - 1) F„_2 = (2.59) 

which follows from the three term recurrence formula (A. 16). 

Note that the five term recurrence formula (2.58) can as well be obtained from the first and second 
order differentiation formulae (2.44) and (2.56) as follows. Take the derivative of equation (2.44), eliminate 
P/.'„ using equation (2.56), and then eliminate the remaining first derivatives Pln+u P'i-n ^'^'^ Pl.n-i using 
equation (2.44) again. 

As the differentiation formulae (2.44) and (2.50), the recurrence formula (2.58) can be also used to get 
step by step Pi-n for all i and n with < i < n, starting from Pi-i (see appendix E). 

2. 5. Recurrence relations with respect to both i and n 

(i) As shown in section 3.3.6 of paper I, SBO polynomials corresponding to different values of i can be 
related together. Thus, P2i+2;2n can be expressed in terms of the polynomials P2i,2e with £ = i, . . . ,n. From 
equations I-(3.58), (2.20) and (2.24), one has for < i < n, 

n-t\n-ij TU+t + e) \m-tj\+i + m 

where, the sum over m can be evaluated as follows. First, the lowest value of m can be taken as I instead 
of sup(i + 1,^), since the term added thereby occurs for t = i only, and it vanishes for m = i. Then, with 
p = n — £, q = m — £, p = £ — i and a=^ + i + £>0, the sum over m reads 



g=0 ^ q=0 

where, using the binomial formula and the beta function [3] 



(2.61) 



Putting all together, one finds for < i < n, 

n-l 



+ P) 



(2.62) 



^2(i+l);2„ = P2i;2n ' 2-"(i +2i){n-i- 1)! ^ JJ^^P2i;2e • (2.63) 

[I I). 



16 



See, e.g., [3] 6.2.1 and 6.2.2. 
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Similarly for the odd polynomials, from equations (2.21) and (2.25), one has for < i < n, 

n-1 2^ 

A(i+l)+l;2n+l = Ai+l;2n+l - 2-"(| + 2i) (n - i - 1)! ^ ^--,^P2i+l;2^+l • (2.64) 

With always the parity relation (2.4), these recurrence relations provide another way of computing Pi-n 
for all < i < n, starting with the boundary values Po;n given by equation (2.26) (see appendix E). Thus, 
for i = 0, one gets for n > 1 [6] 

PiM^) = = 2-3"if2n(\/2a;) - 2-("+i)(n - 1)! ^ H2e{V2x) (2.65) 

2-(2^+i) 

P2;2n+l(x) = P3;2n+l(a;) = 2-(3"+2)/72„+, (^3;) - 2-("+l)3(n - 1)! ^ H2i+l{\f2x) (2.66) 

i.e. expansions in terms of Hermite polynomials of the variable ^/2x, instead of x as given by equations 
(2.20) and (2.21). 

(ii) On the right-hand side of equation (2.63), the sum over £ depends on n through the upper bound n — 1 
of the sum only. Therefore, replacing n by n — 1 in equation (2.63) and eliminating Y^Zi " " " between the 
equation obtained thereby and equation (2.63) yield for < i < n — 1, 

4P2(i+l);2n " 4P2i;2n - 2(n - i - l)P2(i+l);2(n-l) + (2i + 2n - l)P2i;2(n-l) = . (2.67) 

This relation still holds for n = i + 1 since then, it reduces to equation (2.63). From equation (2.64), one 
gets a similar equation for the odd polynomials. Actually, both equations can be written as a remarkable 
single linear four term recurrence relation, with coefRcients independent of x, and which holds for all i and 
n having the same parity and such that < i < n — 2, 

if (- 1)*+" = 1 then 4Pi+2;„ - 4Pi;„ -{n-i- 2)Pi+2-n-2 + (i + n - l)Pi;„_2 = . (2.68) 

For n = i + 2, this relation reduces to a three term recurrence relation. From equations (2.4) and (2.22), one 
finds for i>0, the following expressions which can also be obtained from equations (2.20) and (2.21), 

Pi;i+2 = 2-(^+2) {H,+2 + (1 + 2i)H,) (2.69) 
P,+i;,+3 = P.;.+3 = 2-^'+^^ {H,+3 + (3 + 2i)H,+i) . (2.70) 

More generally, the four term recurrence relation (2.68) with the parity relation (2.4) provide another 
easy way of computing step by step Pi-n for all i and n with < i < n, starting from the boundary values 
of Pi-i given by equation (2.22) (see appendix E). 

2.6. Special value at x = 

The odd polynomials P2i-2n+i = P2i+i;2n+i vanish at a; = 0. For the even polynomials, from equations 
(2.20), (A.ll) and (C.l), one finds for o'< i < n, 

p»-..„(o) = p«„(o) = 2-"t(-ir2-(::')^|f^^^^S)™ 

= (-1)' 2-("-) (i)^(i + 2^)„_^F(-(n - *), i; i + 2^■, 2) (2.71) 
Equation (2.65) has been already given in [6] (3) with different choices of w and W2, see footnote 50. 
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where, F is a terminating hypergeometric function, for which we don't know a closed form in terms of say 
gamma hmction products for all i and n. These values can be computed by recurrence as follows. 

(i) Boundary values for n = i and i = 0. From equations (2.22), (2.26) and (A. 11), one gets 

P2.;2i(0) = 2-''H2M = (-l)^(i), = (-1)^1^ i = 0,l,... (2.72) 

PoMO) = 2-'"if2n(0) = (-l)"2-"(i)^ = (-1)"!^ n = 0, 1, . . . . (2.73) 
Setting for < i < n, 

P2i;2n(0) := (- 1)"2-2("-) (i) . = (-l)"MLp.„ (2.74) 
one finds the boundary values, 

Pi.i = l i = 0,l,... po;n = 2"(i)„ = |^ n = 0,l,.... (2.75) 

(ii) Recurrence relation with respect to n. Setting x = either in the differentiation equation (2.50), or in 
the recurrence formula (2.58), yields for < i < n (setting -f2i;2(i-i)(0) := 0), 

4P2i;2(„+l)(0)+P2i;2„(0)- (n-i)(n + i- i)P2i;2(„-l)(0) =0 (2.76) 

or equivalently from equation (2.74), pi-n satisfies the three term recurrence relation (setting Pi-i-i := 0), 

Pi;n+i - Pi;n - 2(n - i){2n + 2i - l)pi;„_i = . (2.77) 

With the boundary values (2.75), it follows by recurrence that for m > 0, Pi-i+m is a polynomial in i of 
degree [m/2] with positive integer coefficients, e.g.: 

Pi-i+i = 1 pi.i+2 = 8 i + 3 pi-i+s = 24 i + 15 pi-i+i = 192 f + 336 i + 105 . (2.78) 

(ii) Recurrence relation with respect to both i and n. Setting a; = in equation (2.67) yields for < i < n — 1 
(setting P2(i+i);2i(0) := 0), 

4P2(i+l);2„(0) - 4P2i;2„(0) -2{n-i- l)P2(i+l);2(„-l)(0) + {2i + 2n - l)P2i;2(„-l)(0) = (2.79) 

or equivalently with equation (2.74) (setting Pi+i-i := 0), 

2{2i + 1)k+i;„ - Pi-n + 4(2i + l){n-i- l)ft+i;n-i - {2i + 2n- l)pi;„_i = (2.80) 

allowing also to get step by step pt-n for all i and n such that < i < n starting from the boundary values 
(2.75). 

Note that setting a; = in equation (2.44) (resp. (2.56)) allows to compute P/.„(0) (resp. P/.'„(0)), 
knowing Pi-n (0) for several values of n. 

2. 7. Properties of the zeros 

Several properties of the zeros of standard orthogonal polynomials have been known for a long time: e.g., as 
recalled in I-(3.3.7), 

(i) all the zeros are real, simple and located in the support [a, b] of the measure; 

(ii) the zeros of polynomials with consecutive degrees separate each other. 
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Wc saw previously that Hcrmitc SBO polynomials coincide in several special cases with standard or- 
thogonal polynomials, and therefore share their properties, i.e. from equations (2.4), (2.22), (2.26) and 
(2.32), 

P0;0=i?0 = l Pn;n = Pn-l-n ^ H„ n=l,2,... (2.81) 

Po;„(x) =2-5"i^„(^/2a;) Pi.2„+i(a;) = 2-3("+i)if2„+i(V2x) n = 0,l,.... (2.82) 

As a general result proven in section 3.3.7 of paper I: 
(i) for < i < n, Pi-n has at least m distinct real zeros of odd order, with i < m < n; 
(a) Pn-i;m n = 1, 2, . . . has n real and simple zeros. 

It is recalled that the arguments used for standard orthogonal polynomials can no longer be used since 
the scalar product (, )o, defined by equation I-(3.5) (see footnote 12). is not defined by a single integral. 
We observe that all special Hermite SBO polynomials of degree n listed in appendix E have numerically n 
real simple zeros. Nevertheless, we have not been able neither to prove this fact for all n, nor to build a 
counter-example. 

Concerning the relative positions of the zeros of SBO polynomials with a given i and consecutive degrees, 
we only observe numerically that the zeros of the special Hermite SBO polynomials listed in appendix E 
separate each other. Once again, we have not been able to prove a general property. 

3. Laguerre standard block orthogonal polynomials 

The study use the same method and follows the same steps as for the Hermite polynomials in section 2. The 
(generalized) Laguerre polynomials n = 0,1, . . ., see appendix C, correspond to the measure defined by 

a:=0 b:=oo w:= e-'^ x" a > -1 . (3.1) 



3.1. Metric tensor components for a particular weight function W2 
(i) Basic quantity jj^k- With the same interval [a,b], let us choose, 

W2 = w^"'"^ := e"'*^ x" a > -1 > . (3.2) 

This choice is motivated by the application considered in [7] for the special value /i = 2. Our study will be 
also restricted in (ii) below to this special case, since it allows to push the calculation to its end analytically. 
Then, with the same steps as for equation (2.6), using the equations (B.6) and (B.3), one gets 

T(a + l+j)T(a + l + k) 
^jAo^^I^) ■= [^j ' = — -— ^1 -1^ ^"+'>S{j,k,a+l;-iJ, \-iJ. (3.3) 

where, the function S is defined in equation (C.8). Alternately, the result above can be obtained using the 
definite integral given in [5, 10] Another way to get the basic quantity jj^kict, ij) is to compute the scalar 

product {Fs"^ , P/"'')2 where, Fz"\x) is the generating function of Laguerre polynomials given by equation 
(B.17). Then, jj^k{oi,^) is the coefficient of sH'' in the expansion of this scalar product, i.e., 

oo 

(pi") , f/"))^ = ^ ^jA^,f^)sH>' = r{a + 1)(^ + (^ - 2)st - (^ - l)(s + f))-"-^ . (3.4) 

j,k=0 

For /X = 1, i.e. for W2 = w, using equation (C.14), one recovers the orthogonality relation of Laguerre 
polynomials, 7j,fe(a,l) = hjSj^k, with hj given by equation (B.5). This result follows also from equation 
(3.4), using equation (C.l), 

j=0 ^ ^ ^ 3=0 

1^ See, [5] 7.414 4. or [10] 2.19.14.6. p 477, using also the linear transformation formula [9] 15.3.6. 
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(ii) The special case W2 = w^'"^ = e ^^cc". As for Hermite polynomials, when jj, = 2 one gets simple 
expressions. From equations (3.3) and (C.15) setting z = —5, one finds 

7,,,(a,2) = 2-(«+^+^-+'=)^>±li^ (3.6) 

i.e., once again, the same structure with a factor T{c+ j + k) multiplied by some factors depending separately 
on j and k. Alternately, the expression above follows also from equation (3.4), 

(^p(a) ^ p,(a)^^ = r(a + 1)(2 - s - t)-«-i 

= r(a+l)^ ^ (-!)"["" "^V"')2-°-i-"s'"t"-'". (3.7) 

n=0 m=0 ^ / \ / 



Then, setting j ^ m, k = n ~ m and using equation (C.l), one recovers equation (3.6). 
From now on, except explicitly noted, only this special case fi = 2 will be considered. 

Then, the corresponding polynomials Pi-n, < i < n are called the Laguerre SBO polynomials. The special 
case i= \ has already been considered in [7] with different choices of weight functions . (See footnote 52 
in appendix F for interrelations with our polynomial Pi-n-) 

3.2. Connection coefficients between block orthogonal and classical polynomials 

As for Hermite polynomials, the relevant quantities are given in terms of determinants which can be evaluated 
using lemma D. One finds from equation I-(3.14), for < i < n, 

and therefore, from equations I-(3.15), I-(3.7), I-(3.16) and I-(3.17) with &„ given by equation (B.2) [7] 

Hi,^ = 2-(«+i+2") („ _ j)! r(a + 1 + i + n) (3.9) 

hn = ±{-ir2--ml (^^-A Ha + l + i + n) (^^) 
' ' \m-ij T(a+l + i + m) ^ ^ 

4«) = y(_i)™2-"l('"-*')4^±l±^P™ (3.11) 
' n\\m-i}r(a + l+i + m) ^ ^ 

where, according to equation I-(3.13) and (B.2), for i>0, 

P,., (-l)^j!Lf). (3.12) 

The orthogonality and normalization of these polynomials with respect to the scalar product ( , )2 can be 
checked computing Jq°° dxe~^^ x" Pi;n{x) Pi;n'{x) with the expansion (3.10) and equations (3.6) and (C.14), 
thereby recovering equation (3.9). 

See [7] section 2 where, to define the subspace V^, instead of considering the orthogonal polyno- 
mials Qn with respect to (, ), the author computes numerically the "seed" polynomials Q„ := a;" — 
{x^ , a:")/(x° , x^), n = 1,2, . . ., which readily fulfil the constraint {x^ , Qn) = 0. This can be done eas- 
ily because there is only one constraint, i.e. i = 1. 

Equations (3.9) (resp. (3.10)) generalizes the formula given in [7] (7) (resp. (A.l)) for the special case 
i = 1 with diff'erent choices of w and W2, see equations (F.16) and footnote 52. 
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Equations (3.10) and (3.11) should be compared to equations (B.6) and (B.9), respectively. More 
precisely, since w = uig^'"^ following arguments similar to those given for equation (2.30) it follows from 
equations (B.3) and (3.6) that, 

{x'+^ , a;*+'=)(i'") = {x^ , a:'=)(i'"+2*) = T{a + l + 2i + j + k) 

= 2«+i+2* (2^- {i + j)\ L<^] , 2^= (i + fc)! L^l) (3.13) 

where, ( , )('''") denotes the scalar products defined with the weigh function w^'^'°'\ given by equation (3.2). 
Thus, the metric tensor components of the scalar product ( , )(i'"+2») in the basis {x", n = 0, 1, . . .} and 
the metric tensor components of the scalar product ( , )(2.") in the basis {2"(z + n)! n = 0, 1, . . .} are 

equal apart from a common multiplicative factor. Hence, while G-SO with respect to (, )(i'"+2«) of the 
basis {a;", n = 0, 1, . . .} defines the Laguerrc polynomial G-SO with respect to ( , )(2'") of the basis 

{2"(i + n)! n = 0, 1, . . .} defines the polynomial Pi-i+n- In other words, for given i > 0, the linear 

mapping defined by x^ i-^ 2^{i + n)\ i-"'„ , n = 0, 1. ... is such that lI""*"^*^ oc Pi-i+n- Therefore, apart 
from a multiplicative factor c. Pis+n reads as a linear combination of 2"^{i + m) \ L^f^^, m = 0, 1, . . . with 
the same coefficients as in the expansion of L^"^^*-* in terms of the monomials, i.e. from equation (B.6), 

p,„„=.5: ^ ptr^^^r.'r' d^-c^+noiLi:!.. (3.i4) 

n! Tia + 2i + l+ m)\ml 

m=0 ^ / \ . 



With Pi-i+n a monic polynomial, computing the coeHicient of a;" fixes c = (— 1)'2 "n!. Changing i + m into 
m and i + n into n, with now < i < n, one recovers equation (3.10). 

For i = 0, Po;n has to be the unique standard monic polynomial Q2;n corresponding to the weight 
function W2- Since, 'W^''^\x) = 2~°' w^^'°'\2x), with the normalization (B.2), one gets for n > 0, 

Po-,n{x) = 2-"i(f)(2a;) = (-l)"2-"n! L^f )(2a;) . (3.15) 

Thereby, setting z = in equations (3.10) and (3.11) yields the remarkable linear relations between the two 
basic sets of polynomials {Lm \ m = 0, . . . ,n} and {Lm\2x), m = 0, . . . , n} for n > 0, 

Ll^\2x) = i-ir jZ (-ir2- ("+^ + '^|— lW(^) (3.16) 

m=0 

These relations can be checked using the integral already quoted in footnote 18 and the orthogonality relation 
of Laguerre polynomials with hn given by equation (B.5), e.g., one has 

/" dx e-- x'^ Li^H'^^) Lt\x) = (-1)"+™ 2- (« + ^ + "^)»-'" f,^ . (g.ig) 
Jo (n — m)! 

It follows from equations I-(3.10), I-(3.11), (3.10) and (B.7) that, 

Ri;n = -\{i{ct + i) + n{a + n)) (3.19) 

Si-^ = ^(^i^[2a + i) + (2n^ + 2(a - 3)n - 3(a - 1)) i{a + i) + n{n - l){a + n){a + n - 1)J . (3.20) 

What is new in the Laguerre case with respect to the Hermite one is that a shift of i in the powers of x 
can be interpreted as a shift in the parameter a. This is why in the Hermite case, this argumentation can 
be used with i = only. 
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For general properties of SBO polynomials about the projection operators, the integral representations 
and the normalizations, see sections 3.3.1, 3.3.2 and 3.3.4 of paper I, respectively. 

3.3. Differentiation formulae 

For the scalar product (, ), the weight function w satisfies the differential equation (B.ll), which implies 
the relations (B.12) and (B.13). Similarly, for the scalar product (, )2 corresponding to the weight function 
yj(/'>a) (defined by equation (3.2), one has for any polynomials / and g, 

= -a- 1)4'''"^ (3.21) 

{f,xg')2 = - {xf -{i^x-a-l)f,g)^ (3.22) 
{xf" -{fix-a-l)f',g)^ = -{f',xg')^ = {f, xg" -{ixx-a- l)g') ^ . (3.23) 

(i) For < ^ < i — 1, the polynomial xL^' of degree ^ is in the subspace Vi. Therefore, xL'}' and Pi-n for 
i <n are orthogonal with respect to ( , ). Prom equation (B.12), one has for < ^ < i — 1 and < i < n, 

Hence, the particular polynomial xP-.^—xPi-n of degree n+1 is orthogonal to the subspace Vi with respect to 

( , ), and thus reads oa-.njPi-.j- Computing the coefficient of a;""*"^ in a;P/.„ — xPi-n fixes ai-n^n+i ~ ^1- 

Now, with respect to the scalar product ( , )2, one has from equations (1.2) and (3.22), for i < j < n + 1, 

= -{xPlj - (m - l)xPi,j , Pi,n)2 - (a + ^)Hi;nhn ■ (3-25) 

In the special case fi = 2, \t turns out that xPlj — (/x — l)xPi-j is the particular polynomial xPlj — xPi-j, 

introduced above for j = n (i.e. of degree j + 1 in the subspace V^), and which reads J2k^i oti-j^kPi-k- Thus, 
for = 2, considered henceforth, one has 

ai;n,jHi.j = - (ai-j^n ^ h,n + (a + l)'^j,n) Hi-n ■ (3.26) 

k=i 

It follows that for j = i, . . . ,n — 2, ai-n.j vanishes. It remains to consider the cases j = n — 1 and j = n. 
With ai-n-i,n = —1 and equation (3.9), one gets 

H 1 1 

Ki;n{a) := a»;n,n-i = ^ = -{n - i){a + i + n) = --(a + 1) • (3.27) 

Hi;n-1 ^ ^ 

Finally, for /* = 2 and < i < n, one finds the first order differentiation formula [7] 

XP;.^ = + (x - i(a + 1))P,;„ + K^■,nia)Pr,n-l (3.28) 

which still holds for n = i, since then Ki-n(oi) vanishes. This equation has to be compared with the differen- 
tiation formula (B.16). For i = 0, using equation (3.15), one gets 

2x4")' = (n + l)4"^^i + ix-a- 1)4") - (n + a)4"_\ (3.29) 

This generalizes the formula given in [7] (3) for the special case i = 1 with different choices of w and 
W2, see footnote 52. 
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which follows from, e.g., adding the two differentiation formulae (B.14) and (B.16). 

With the boundary value of Pi-i given in terms of Laguerre polynomials by equation (3.12), the differ- 
entiation formula (3.28) allows to get step by step Pi-n for all i and n with < i < n. A table of these 
polynomials for z = 0, 1, 2 and n = i,i + 1, ... ,i + A is given in appendix F. In particular, setting n = i + 1 
and using equations (3.12), (B.2) and (B.16), one finds for n = i > 0, 



P;^+l = + -{a + l + 2^) = (-l)'+^(.: + 1)! l'^, + (-1)' - z!(a + 1 + 2z) L^") . (3.30) 



This result can be obtained directly setting n = ?' + 1 in equation (3.10). 

(ii) Following the same method again, for < ^ < i < n, the polynomial xL"" — {x — a — l)i^' of degree £ 
is orthogonal to Pi-n with respect to the scalar product ( , ), and from equation (B.13), one has 







, xLf" -{x-a- 1)4")') = {xP^^ -{x-a- l)?^^^ , L^) . (3.31) 



Hence, the particular polynomial xP"^ — {x — a— l)f'/.„ of degree n belongs to V^- and reads 7i;n,jPi;j- 
Computing the coefficient of in .tP/.'„ — {x — a — l)-P/.„ fixes 7i;„,„ = — n. With the scalar product ( , )2, 
one has from equations (1.2), (3.22) and (3.23), for i < j < n, 

7i;n,jHi.j = [Pi-J , xP^.n — {x — a — l)Pj.„)2 

= {Pi;j ,xPi[n-{t^x-a- 1)^„)2 + (/X - l){h,i , xPl^n)2 

= (xPl'.j - ((m -l)x-a- 1)PQ , Pr,n)^ - l^ixPlj -ill- l)xPi.,j) , Pi;„)^ 

-{|I-l){a + l)H^■,n6j,n. (3.32) 

In the special case ^ = 2, once again, xP^'j — ((/i — l)a; — a — l)P/.j (resp. xP^.j — {fj, — l)xPi.j) is the 
particular polynomial xP,'!^ — {x — a — ^)P'' j (rcsp. xPlj — xPi-^j), introduced above for j = n (i.e. of degree 
j (resp. j + 1) in the subspace Pj^), and which reads J2k=i li;j,kPi;k (resp. J2kt^i '^i;j,kPi;k)- Thus, one has 

Ji;n,jHi-j = f 7i;j,n ^ Sk,n - "^Cti-J.n ^ 5k,n - {a + l)Sj,n] Hi-n . (3.33) 
k=i k=i 

For n = i, and therefore j = i, it is already known that -ji-i^i = —i. For n > z + 1, it follows from above 
that for n > i + 2 and j = i, . . . ,n — 2, Ji;n,j vanishes. It remains to consider the case j = n — 1. From 
C(i-n-i,n = —1 and equation (3.27), one gets 

7i;n,n-l = -2Q;»;n-l,n gj. = 2Ki;„{a) . (3.34) 

Finally, in the case = 2, one finds the second order differentiation formula for < i < n [7] 

xP!!„ + {a + l- x)Pl^^ + nPi,n = 2Ki,n{a)Pi,n-i (3.35) 

where, Ki-^n{oi) is defined by equation (3.27). For n = i, the relation still holds, since then Ki;n(Q;) vanishes, 
and using equation (3.12), the result coincides with the differential equation (B.15). For i = 0, using equation 
(3.15), the relation follows, e.g., from equations (B.14) and (B.15). 

3.4- Recurrence relation with respect to the degree n 

This generalizes the formula given in [7] (6) for the special case i = 1 with different choices of w and 
w^, see footnote 52. 
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Taking the derivative of the first order differentiation formula (3.28), and eliminating a;P/.'„ with the second 
order differentiation formula (3.35) yield 

-Pi;n+i + ^(a - + (n + + \in ' i){a + l + i + n)iP(.^_, - 2Pi;„_i) = . (3.36) 

Then, multiplying both sides of this equation by x and using again equation (3.28) to eliminate the first 
order derivatives Pl-n+n ^tn ^''and -P/.„_i give the five term linear recurrence formula for < i < n, 

Pi;n+2 -{x- l)Pi;n+i + (^{a + 1 + 2n)x - 2Ki;„(a) - i(a^ + a + 2n)^Pi.n 

- Ki;n{0i){x + l)Pi;n-l - Ki;n{a)Ki;n-l{a)Pi;n-2 = (3.37) 

where, Ki-n{ct) is defined by equation (3.27). A direct calculation shows that this relation still holds for 
n = i + 1 and n = i, where it reduces to a four and three term recurrence formula, respectively. For n = i, 
using equation (3.30), one finds for i>0, 

Pi,i+2 = ix- + ^(a^-a-2- 2i)L\"^ . (3.38) 

As the differentiation formula (3.28), the recurrence formula (3.37) can be used to get step by step Pi-n 
for all i and n with < i < n starting from Pi-i and Pi;i+i given by equations (3.12) and (3.30), respectively 
(see appendix F). 

3. 5. Recurrence relations with respect to both i and n 

(i) Setting j = i + 1 in equation I-(3.58) and using equations (3.10) and (3.11), one can express Pi+i-n in 
terms of the polynomials Pi-i with i < £ < n. The calculation is very similar to what was done for Hermite 
polynomials in section 2.5. One finds 



- /„_,x r(a + 2 + i + n) " /^_A ^^j 

and using equations (2.61) and (2.62) (with p = n — i, p = £ — i and cr = a + l + i + f>0)to compute the 
sum over m yield for < i < n, 

n-l I 

Pi+i;n = hn " 2-"(a + 1 + 2i){n -i-l)\Y, 773^^^;^ • (3-40) 

This recurrence relation provide another way of computing step by step Pi-n for all i and n with < i < n, 
starting with the boundary values Po;n given by equation (3.15) (see appendix F). Thus, for i = 0, one gets 
for n > 1, 

Pi;n = (-l)"2-"L««)(2a;) - 2-"(a + l)(n - 1)! ^(-l)^4"^(2a;) (3.41) 

i.e. an expansion in terms of Laguerre polynomials of the variable 2x, instead of x as given by equation 
(3.10). 

(ii) On the right-hand side of equation (3.40), the sum over £ depends on n through the upper bound n — 1 
of the sum only. Therefore, replacing n by n — 1 in this equation and eliminating J2e=i ' ' ' between the 
equation obtained thereby and equation (3.40) yield for < z < n — 1, 

2Pi+i;„ - 2Pi;„ -{n-i- l)Pi+i;„-i + {a + i + n)Pi;„_i = . (3.42) 

Keeping the term P!.^ provides another first order diflFerentiation formula. 
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This relation still holds for n = where it reduces to a three term recurrence recovering equation (3.40). 

This four term linear recurrence relation provides an easy way of computing step by step Pj:„ for all i 
and n with < i < n, starting from the boundary values of P^i given by equation (3.12) (see appendix F). 

3. 6. Special value at x = 

From equations (3.10), (B.8) and (C.l), one finds for < i < n, 

Pi;n{0) = (-l)'2-("-'5 (a + l)i{a + 1 + 2i)„_iF(-(n - i),a + 1 + i;a + 1 + 2i;2) (3.43) 

where, F is a terminating hypergeometric function, for which we don't know an explicit expression for all i 
and n. These values can be computed by recurrence as follows. 

(i) Boundary values for n = i and i = 0. From equations (3.12), (3.15) and (B.8), one gets 

P,,(0) = (-l)'z!if\0) = (-l)*(a + l), ^ = 0,l,... (3.44) 

Po;„(0) = (-l)"2"n!4")(0) = (-l)"2"(a+l)„ n = 0,l,.... (3.45) 

Setting for < i < n 

Pi;40) := (-1)" 2-("-^) {a + l)ipg (3.46) 

one finds the boundary values, 

pg = l i = 0,l,... pg = (a + l)n n = 0,l,.... (3.47) 

(ii) Recurrence relation with respect to n. Setting a; = in the differentiation equation (3.28) yields for 
< i < n (setting Pi-i-i{0) := 0), 

4P,;„+i(0) + 2{a + l)Pi;„(0) - (n - i){a + i + n)P,;„_i(0) = (3.48) 

or equivalently from equation (3.46), p\"J^ satisfies the three term recurrence relation (setting p\"^_i := 0), 

pS+1 - (" + -in-z){a + i + n)p'^_, = . (3.49) 

With the boundary values (3.47), it follows by recurrence that for m > 0, Pi"+„i is a monic polynomial in a 
of degree m with positive integer coefficients, and also a polynomial in i of degree [m/2], e.g.: 

f Sli = " + 1 P&^2 = a' + 3 a + 2(1 + i) pfg^, = (a + 1) (a^ + 5 a + 6(1 + i)) 

p\"^^^ = a^ + lOa^ + (35 + 12 i) o? + 2(25 + 18 z) a + 12(1 + i)(2 + i) . (3.50) 

(Hi) Recurrence relation with respect to both i and n. Setting x = in equation (3.42) yields for < i < n— 1 
(setting P.+i;.(0) := 0), 

2Pi+i;„(0) - 2Pi;„(0) - (n - i - l)Pi+i;„_i(0) + {u + i + n)Pi;„_i(0) = (3.51) 

or equivalently with equation (3.46), 

2(2a + 1 + z)p^\.^ - p^J +2{a + l + z){n-t- l)plXl,n-i + * + n)p^^_, = (3.52) 

2^ See equation (3.56) below, pi-n defined for Hermite polynomials by equation (2.74) and p\"^ are related 
bypi;„ = 2"->y\ 
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allowing also to get step by step p^^^ for all i and n such that < i < n starting from the boundary values 
(3.47). 

Note that setting a; = in equation (3.35) (resp. its derivative with respect to x) allows to compute 
P/.„(0) (resp. P/'„(0)), knowing Pi;„(0) for several values of n. 

3. 7. Properties of the zeros 

From equations (3.12) and (3.15), the particular Laguerre SBO polynomials, 

P„;„ = £i«) Po;„(a;) = 2-"Li«)(2x) n = 0,l,... (3.53) 

share the properties, especially of the zeros, of the standard orthogonal polynomials. 
It has been shown in section 3.3.7 of paper I that: 

- for < z < n, Pi-n has at least m distinct non-negative real zeros of odd order, with i <m < n; 

- for n> 1, Pn-i-n has n non-negative real and simple zeros. 

We observe that all the special Laguerre SBO polynomials of degree n listed in appendix F have numerically 
n real non-negative simple zeros. Nevertheless, we have not been able neither to prove this fact in general, 
nor to build a counter-example. 

Concerning the relative positions of the zeros of SBO polynomials with a given i and consecutive degrees, 
in the Laguerre case one readily finds a counter-example, e.g., with a = 0, from equations (F.6) and (F.7), 
the zeros read 

P,.2=a;2-^x+i a;i = « 0.219 = « 2.281 (3.54) 

2 2 4 4 

Pi.3 = a;3-5x2 + 5x-l t/i = 2 - ^3 « 0.268 t/2 = 1 2/3 = 2 + Vs « 3.732 (3.55) 

where, xi < yi < y2 < X2 < y3, and therefore the zeros do not separate each other, as for standard orthogonal 
polynomials. Once again, we have not been able to prove a general property. 

3.8. Interrelations between Hermite and Laguerre cases 

The interrelations between Hermite and Laguerre polynomials are recalled in appendix B (item 6). Using 
the relation (B.19) (resp. (B.20)), the expansion (2.20) (resp. (2.21)) for the even (resp. odd) polynomials 
in the Hermite case, and the expansion (3.10) for the Laguerre case, setting a = —5 (resp. 5), one finds for 
0<i<n, 

= = Plir'^'''' ''=--\x^) (3.56) 

^(a^) = pst:'^:u^) =^p^'"' (3.57) 

i.e. the same interrelations as between Hermite and Laguerre polynomials. Let us recall that, while the choice 
of i determines the subspace Vi of polynomials of degree at most i — 1, the complementary and orthogonal 
subspace V^, equation (1.3), depends on the choice of the scalar product (, ). Therefore, ■pj-^'^™^^'^)-^ ^ 

p^Laguerre a- j)^- ^^^^ p^(Laguerre a-2)-L ^^^^ ^ priori, different SUbspaCCS. 

Relations (3.56) and (3.57) allow a cross check of many similar formulae obtained in sections 2 and 3, 
and also of the explicit expressions of the Hermite and the Laguerre SBO polynomials listed in appendices 
E and F, respectively. 

Conclusion 

In a first paper [1] (paper I), we have defined general block orthogonal polynomials. Their properties have 
been investigated, pointing out what can be generalized and ascribed to Euclidean vector space, endowed 
with two distinct scalar products, and to polynomial algebra. A general study of the special case of standard 
block orthogonal (SBO) polynomials, has been made, emphasizing the similarities to and the diff'erences from 
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the study of the standard orthogonal polynomials. For a given « = 0, 1, . . ., these real polynomials of degree 
n, Pi-n, n = + . . ., are orthogonal to the «-dimensional subspace Vi of polynomials of degree less than i 
with respect to a first Euclidean scalar product ( , ) . In addition, they are mutually orthogonal with respect 
to a second Euclidean scalar product (, )2- They form an orthogonal basis with respect to ( , )2 of the 
subspace of codimension i, orthogonal to Pi with respect to ( , ), and complementary to Pi with respect 
to the space of polynomials of any degree. Based on the Gram-Schmidt orthogonalization (G-SO), we have 
determined these polynomials uniquely, apart from a constant multiplicative factor for each polynomial, in 
terms of determinants which entries are the metric tensor components of both scalar products considered. 

In this second paper, we apply the general results obtained in paper I, assuming the first scalar product 
( , ) corresponds to the Hermite (resp. Laguerre) classical polynomials. In both cases, it turn out that for 
a particular choice of the second scalar product ( , )2, we get explicitly the connection coefficients between 
the Hermite (resp. Laguerre) polynomials, and what we call the Hermite (resp. Laguerre) SBO polynomials. 
These polynomials satisfy all the general properties derived in section 3.3 of paper I for the projection 
operators, the integral representations and the normalizations. Furthermore, we establish recurrence relations 
with respect to n and i, and differentiation formulae, allowing a step by step determination of these SBO 
polynomials. As for the classical orthogonal polynomials, these relations are based on differential equations 
satisfied by the weight functions w(x) and W2{x) defining, respectively, the two scalar products. Partial 
properties for the zeros of these polynomials are also obtained. We point out the differences with the 
classical orthogonal polynomials, basically due to the fact that two, instead of one, distinct scalar products 
are involved. As the Hermite and the Laguerre polynomials, the Hermite and the Laguerre SBO polynomials 
can be related together. In the Hermite case, the measures considered being symmetric with respect to the 
origin, additional parity properties occur. 

This paper recovers and improves results already obtained by Giraud et al [6- 8] for i = 1, and partly for 
i = 2. In particular, we get new recurrence and differentiation formulae. We generalize these results making 
a systematic study for all non- negative integer values of i. As recalled in paper I, Pj;„, n = i, i + l, . . . provide 
a basis of the subspace P^ ^ suitable for applications requiring to take into account i linearly independent 
constraints which can be recast into an orthogonality relation with the subspace Pi with respect to (, ). 
(We have shown in paper I that an extension of similar methods to constraints associated with different 
scalar products is not straightforward, and may have no solution, a unique solution or an infinite number of 
solutions.) These constrained polynomials have been first considered in [6-8] with i = 1 io implement the 
matter conservation in view of application of the Hohenberg-Kohn functional in some problems of nuclear 
physics. 

Open questions still remain. Why can we perform analytical calculations only for very special choices 
of the weight function W2 (precisely those made in the applications above), once w is chosen as the classical 
Hermite (resp. Laguerre) weight function? Is there any 'integrability' structure to discover? Nevertheless, 
let us recall that for any proper weight functions, determination of the SBO polynomials is a well-defined 
problem. G-SO is well-suited for numerical calculations and it is not necessary for applications to have 
analytical expressions. We have not been able to get closed formulae for the connection coefficients between 
Pi-n and a;™ (we observe in all the special values of Pj:„ given in appendices E and F, alternate signs in the 
expansion in powers of a;). Can we establish that, both in Hermite and Laguerre cases, Pi-^n has n distinct 
real zeros? Similarly, do the zeros of the Hermite SBO polynomials with a given i and consecutive degrees 
separate each other? (We establish this is not true for the Laguerre case.) 

As a byproduct, for completeness and to underlined the similarities to and the differences from the 
classical polynomials, we recover in appendices well-known properties of Hermite and Laguerre polynomials 
in a unified way (including the expansions of monomials in terms of classical polynomials, not given in the 
standard textbooks). In addition, also reported in appendix, we have been led to consider several terminating 
generalized hypergeometric series with one or more variables often set to unity. 

The Gegenbauer and Jacobi polynomials will be similarly considered in a forthcoming paper. 

We have already underlined in the conclusion of paper I (see footnote 44 of paper I), that the general 
formulae established therein can be used to defined block orthogonal polynomials for any given i-dimensional 
subspace 5,, possibly different from Pi (see footnote 51). Then the block orthogonal polynomials would not 
be the standard ones and require a new similar study. 
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Appendix A. Hermite polynomials 

It is shown in appendix I-A that all the relevant quantities (Z„, /i„, Qn, bm.n, ■ ■ •) can be expressed in 
terms of determinants with the metric tensor components as elements, except possibly for the last raw, 

Qj^k ■■= {x^ , x'') = Cj+k ■■= / dxwix) x^+'' . (A.l) 

J a 

These general results are now applied to the classical orthogonal polynomials, i.e. for particular choices of the 
measure defining the scalar product. Thereby, well-known results are recovered using a unique method based 
on the calculation of some special determinants with gamma functions as elements (see lemma D in appendix 
D) and also the evaluation of some finite sums involving binomial coefficients and ratios of gamma functions 
(sec lemmas C.l and C.2 in appendix C). In addition, since the weight fimctions w considered satisfy some 
differential equations, integrations by part with vanishing boundary conditions yield some relations between 
scalar products {this is still valid for particular SBO polynomial, see sections 2.3 and 3.3.). 

(1) Definition. Hermite polynomials Hn, n = 0, 1, . . ., are defined with the measure and normalization 

a := —DO b := oo w := e~^^ (^-2) 
fc„ := 2" n = 0, 1, . . . . (^.3) 
The measure being symmetric with respect to the origin, Hn{—x) = (— l)"i/„(x). 

(2) Connection coefficients. Prom equation (A.l) and the Euler integral the nonzero moments read 



C2j 



/OO 
dxe--\^^ =r{^+j) j = 0,l,.... (A4) 
-OO 



Therefore, all the determinants which occur in equations I-(A.17)-I-(A.28) can be obtained with lemma D 
and the duplication formula (see footnote 6). For n = 0, 1, . . . and m = 0, . . . , n, one finds 

n n 

Zl^^ = X{3\V{\+3) Z(°) = nj!r(|+i) (A5) 
/i„ = 2"n!7r5 (A6) 

a2™+i,2„+i = fc2.+i(-l)"-"' 3^^U ") = (-l)"-'"2^"-+^ (^.8) 

r(|+m)VW (2m + 1)! (n - m)! 

[»/2] , 

Hn=y {-IT ^, (2a;)""^'" (A.9) 

m=0 ^ ' 

r„=0 Sn = --n{n-\) (^.10) 



.(2n)! 

See, e.g., [4] 10.13 (1) and (4). 



H2n{0) = (-l)"22"(i)„ = (-1)"^ H2n+l{0) =0. (^.11) 



26 

27 See, e.g., [3] 6.1.1. 

2^ See, e.g., [4] 10.13 (4) for hn, (9) for the expansion (A. 9), (4) for Vn only (see footnote 47 in paper I), 



(15) for Hn{0) and (5) for the three term recurrence formula. 
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Similarly, the connection cocfRcient bm,n, not given in the standard textbook already quoted, can be obtained 
from equations I-(A.IO), I-(A.22), I-(A.28), (D.2), (C.3) and the duplication formula (see footnote 6). For 
m = 0, . . . , n one has 

-2™r(i + n) /n\ ^_2„ (2n)! 



= 2-^"— ^ =2-^"- (A12) 

r(i+m)V™/ (2m)!(n-m)! ^ ^ 

. _ o-(2,n+i) r(| + n) fn\_ (2n + 1)! 

02m+l,2n+l — ^ ' /o , T — ^ 77^ , .x, . 77 [A.L6) 

r(|+m)VW (2m + l)!(n-m)! 

ml in — 2m)! 

m=0 ^ ' 

These expressions of bm,n can also be obtained directly with the following lemma, setting pk = (— 2^*^ r(^ + 
k) and = 

Lemma A. With pj and aj some non-zero complex numbers, the matrices A and B defined by, Aj^k '■= 
{—iyQ)pk/o'j Bind Bj^k '■= {~^y {'j)(^k/Pj for j,k = 0,1 .. ., are inverse one from the other. 

Proof. Since the binomial coefficient (^) vanishes if fc < j, the matrices A and B are triangular with zeros 
under the diagonal and nonzero diagonal elements. Then, this known lemma is a direct consequence of 
the binomial formula. Indeed, AB is also a triangular matrix and one has 

(3) Three term recurrence formula. From equations I-(3.55), (A. 3) and (A. 10), one has for n = 0, 1, . . . 

Hn+l = {AnX + Bn)Hn-CnHr,-l A„ = 2 B„ = C„ = 2n . (A16) 

(4) Differentiation formulae. With w defined by equation (A. 2) and / and g some polynomials, one has 

w' = -2xw (^-17) 

{wf)dg = - d{wf)g = -{f' -2xf,g) (A18) 

-OO J — OO 

(/" - 2xf' , g) = -if , g') = if , g" - 2xg') . (A.19) 

These equations can be used in two ways to generate differentiation formulae for n = 0, 1, 

• First method 

(i) As any polynomial of degree n—1, H'^ = J2"Zo '^j^j where, from the orthogonality relation I-(A.2) and 
equation (A. 18), 

ajhj = {Hj ,K) = - {H'. - 2xHj , if„) . (A.20) 

2® See, [3] table 22.12 where, the numerical values of bm,n, and also am,n, are listed up to n = 12, but 
the general expressions (A.12)-(A.14) are not given there. Formulae (A. 12) and (A. 13) can be found in [11] 
chapter 6, exercise 6 p 339. 

See, e.g., [16] p 3, going from equation (1.8) to equation (1.91). 

31 See, e.g., [4] 10.13 (5). 
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Since H'j — 2xHj is a polynomial of degree j + 1, the sealar product above, and thus aj, vanishes for 
j = 0, . . . ,n — 2. Computing the coefficient of x"~^ in H!^ fixes a„_i = nkn/kn-i- One gets 

H'^ = 2nff„_i . (A21) 

Following the same steps, first order differentiation formula for x^H^, with p = 0, 1, . . ., can thereby be 
obtained. Thus, for p = 1, considering the polynomial xH!^ of degree n, one finds 

xH!^ = nHn + 2n{n - l)F„-2 • (^.22) 

(ii) Similarly, being a polynomial of degree n, H'^ — 2xH^ = ^2^=0 ^j^j where, from equation (A. 19), 

ajhj = {Hj , Hll - 2xH^) = {H'^ - 2xH'. , H„) . (^.23) 

Since H'- — 2xH'j is a polynomial of degree j, aj vanishes for j = 0, . . . , n — 1. Computing the coefficient of 
a;" in H'^ — 2xH'^ fixes a„. Thereby, iJ„ satisfies the second order differential equation (2.36) 

A{t)H'^ + B{x)H'^ + XnHn = A{x) = 1 B{x) = -2x A„ = 2n . (A.24) 

• Second method (permutating n and j in the relations above) 

Only this method can be extended to the SBO polynomials as shown in section 2.3. 

For j = 0, . . . ,n, the polynomial H'j of degree j — 1 is orthogonal to and with equation (A. 18), 

= {H„ , i/j) = « - 2xH„ , H,) . (A.25) 

Hence, the particular polynomial i/^ — 2xHn of degree n + 1 is orthogonal to the subspace Vn+i, and thus 
it reads a„+ii/„+i. Computing the coefficient of in iJ^ — 2xHn fixes a„+i. One gets 

K - 2xHr, = -Hr,+i . {A.26) 

Following the same steps, first order differentiation formula for x^H^, with p = 0, 1, . . ., can thereby be 
obtained. Thus, with p= 1, considering the polynomial xH'j of degree j, one finds 

xir; = -i7J„+2 + (2^2-1 -n)if„. (A27) 

Actually, all these differentiation formulae and the three term recurrence formula are not independent. 
Using simple algebra (substitutions and derivations), it can be checked that: 

(i) all these differentiation formulae follow from the first order differentiation formula (A. 21) and the three 

term rccmrcnce formula (A. 16); 

(ii) more generally, among the equations (A. 16), (A. 21), (A. 24) and (A.26), any pair of equations imphes the 
other two equations. In particular, the first order differentiation formula, e.g., (A.26), and the second order 

differential equation (A. 24) yield the recurrence equation (A. 16). This kind of method is used in section 2.4 
to find a five term recurrence formula with respect to the degree for the Hermite SBO polynomials in the 
special case W2 ■= e~^^ ; 

(iii) note that the differentiation formulae (A. 22) and (A. 27) relate Hermite polynomials having the same 
parity. From both equations, one gets the recurrence formula, 

Hn+2 = 2(2a;2 - 2n - l)i?„ - 4n(n - l)i?„_2 (A.28) 

which can be obtained also iterating the three term recurrence formula (A. 16). 

32 See, e.g., [4] 10.13 (14). 

33 See, e.g., [4] 10.13 (12). 
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(5) Generating function. The most common generating function reads 

oo _. 

F,{x) := ^ - Hn{x) = e^^'-'" . {A.29) 



n=0 



Appendix B. Laguerre polynomials 

(1) Definition. (Generalized) Laguerre polynomials n = 0,1,... are defined with the measure and 
normalization 

a:=0 6:=oo w := e'^o;" a > -1 (B.l) 

k„:=^-^ n = 0,l,.... {B.2) 
n! 

(2) Connection coefficients. Prom equations (A.l) and the Euler integral (see footnote 27), the moment reads 

POO 

Cj= dxe-''x"+^ =T{a + l+j) j = 0, 1, . . . . (B.3) 

Then, all the determinants which occur in the general formulae I-(A.5)-I-(A.10) can be obtained with lemma 
D. For n = 0, 1, . . . and m = 0, . . . , n, one finds 

n 

Z„ = ]lj'.r(a+l+j) (S.4) 

n! 

(-1)" r(a + l + n) fn\ j-l)^ f a + n\ 
n! L [a + 1 + m) \m J m! \n ~ m J 

rn = -n{a + n) s„ = ^ n(n - 1)(q + n)(a + n - 1) {B.7) 

Similarly, the connection coefficient 6m, «, not given in the standard textbook already quoted, can be obtained 
from equations I-(A.IO), (D.2) and (C.3). For m = 0, . . . ,n one has [3] 



™! ^ V^y r(a + l + . 

^ ^ T{a + l+m)\mJ ^ ' \n-mj ^ ' 

34 See, e.g., [4] 10.13 (19). 

3^ See, e.g., [4] 10.12 (1) and (2). According to [3] 22.2.12, a is written between parentheses to avoid any 
confusion with a power index. 
3^ See, e.g., [4] 10.12 (2) for /i„, (7) for am,n, (2) for r„ only (see footnote 47 in paper I) and (13) for 

3^ See, [3] table 22.10 where, the values of bm,n, and also am,n, are listed only up to n = 12, but the general 
expression (B.9) is not given. 
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This value of bm.n 

can also be obtained directly with lemma A, setting pj = T{a + 1 + and aj = 

r(a + l+ j). 

(3) Three term recurrence formula. From equations I-(3.55), (B.2) and (B.7), one has for n = 0, 1, . . . 

(n+l)A„ = -l {n+l)Bn = a + 2n+l (n + l)C„ = a + n. (B.IO) 

(4) Differentiation formulae. With w defined by equation (B.l) and / and g some polynomials, one has 

{xw)' = -{x~ a~l)w (B.ll) 

if , xg') = xwfdg = - d{xwf) g = - {xf -{x-a-l)f,g) (B.12) 
Jo Jo 

[xf" -{x-a- 1)/' , g) = -if , xg') = -{xf , g') = (/ , xg" -{x-a- l)g') . (S.13) 

As for the Hermite polynomials in appendix A 4-, these relations can be used in two ways to generate 
differentiation formulae for n = 0, 1, 

• First method 

(i) Considering the polynomial xLI"'' with equation (B.12), one finds 

xLi"^' = nL^-^ - (a + n)4"_\ . (B.14) 

(ii) Considering the polynomial xL^n^" — {x — a — with equation (B.f3), one gets the second order 
differential equation (2.36) 

A(x)L(f)" + B(x)L(^")' + A„i(f) =0 A{x)=x B{x) =a + l-x A„=n. (B.lb) 

• Second method (permutating n and j in the relations above). Only this method can be extended to the 
SBO polynomials as shown in section 3.3. 

For j = 0, . . . , n — 1, considering the polynomial xL^"^', one finds 

xLl^^' = {n + 1)L^^^, + (x - Q - 1 - n)4") . (B.16) 

As for the Hermite polynomials, it can be checked that: 

(i) all these differentiation formulae follow from the first order differentiation formula (B.14) and the three 
term recurrence formula (B.IO); 

(ii) more generally, among the equations (B.IO) and (B.14)-(B.16), any pair of equations implies both other 

equations. In particular, the first order differentiation formula (B.16) and the second order differential 
equation (B.15) yield the recurrence equation (B.IO). This kind of method is used in section 3.4 to find a 
five term recurrence formula with respect to the degree for the Laguerre SBO polynomials in the special case 
W2 := e~'^^x°'. 

(5) Generating function. The most common generating function reads 

oo 

i^W(x):=^4")(a;)2" = (l-z)-"-^exp^^ \z\ < 1 . {B.17) 

n=0 ^ ^ 
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See, e.g. 


[4] 10.12 (3). 
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See, e.g. 


[4] 10.12 (12) 


40 


See, e.g. 


[4] 10.12 (10) 
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See, e.g. 


[4] 10.12 (12) 


42 


See, e.g. 


[4] 10.12 (17) 
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(6) InterrelaMons between Hermite and Laguerre polynomials. For any polynomials / and g, a change of 
variable in the integral yields 

r dxe-^'fix)g{x)= r dxe-''x-^f{^)g{^)^ H dxe'^x^ ^^^^ . {B.18) 

J-oo Jo Jo VX y/x 

It follows that forn = 0, 1, . . . 

Li-^\x) = H2n{V^) H2n{x) = (- l)"22»n! 4" ^^x^) (B.19) 

Lt^\x) = ^^^^ H2n+i{x) = {-ir2'-+'n\xLi^\x'). {B.20) 
y/x 

Appendix C. Pochammer symbol and some related relations 

Definition and properties of the Pochammer symbol or rising factorial {z)n are detailed in [13, 14] With 
n a non- negative integer and z some complex number, one has , 

■= ^%y^ = {~n) =z{z + l)---{z + n~l) = i-in-z - n + 1)„ . (C.l) 

The Pochammer symbol satisfies the binomial formula, i.e. with n a non-negative integer and z and w some 
complex numbers , 

(^ + «;)„ = f^(")(^),H„_,. (C.2) 

i=o 

Lemma C.l. With j and k non-negative integers and c some complex number not zero or a negative integer, 
one has 

where, the binomial coefEcient (^) vanishes if j < k. 

Proof. Setting n = k, z = c + j and w = —c — k + 1 in the binomial formula (C.2) yields 



£=0 

Then, it follows from equation (C.l) that. 



ij-k + l)fe = J2 (^) (c + Jh {-c-k + l)k-e ■ 



(C.4) 



thereby recovering equation (C.3) provided c + j ^ 0, —1, This identity can be proved also directly by 

recurrence on k. 



43 See, e.g., [4] 10.13 (2) and (3). 

See, e.g., [13] appendix I or [14] section 2. 

45 See, e.g., [3] 6.1.5, 6.1.21 and 6.1.22. 

46 See, e.g., [14] section 2.6. 
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Lemma C.2. With c, z and w some complex numbers, c not zero or a negative integer, and j and k 
non-negative integers, let, 

minO.fe) , , , ^ 

S,{j,k,c;z,w):= ^ (^^j y (,^). (1 + (1 + (C.6) 



p=0 

^ ^ VpA?; r(c + p) r(c + 



S{J, k, c; z, H := F ^-j, -k; c; ) (C.S) 

where, F is an hypergeometric function. Then, one has 

S{j,k,c;z,w) ^ S{k,j,c;w,z) = Si{j,k,c;z,w) = S2{j,k,c;z,w) . (C.9) 



Proof. From equations (C.l), in terms of the Pochammer symbol, one has 

r(c) ^ {c)p p\ \{i + z){i + w)j 

= S{j,k,c;z,w) (C.IO) 

r(c) (c)p(c)g p! g! 

F2(c, -j, -fc,c,c;-2;,-u;) (<^-ll) 



r(c) 



where, F2 is a special case of a generalized hypergeometric series of two variables listed by Horn and reducible 
to an hypergeometric function [12] 



F2(c, -j, -k, c, c; -z, -w) = (1 + (1 + F -j, -fc; c; ^ , ^ . . (C.12) 

V {l + z){l + w)J 

The symmetry property (j, -z) is obvious. This completes the proof of equation (C.9). 

Corollary C. Special values of S{j, k, c; z, w), see lemma C.2 equations (C.6)-(C.9): 

^(.,fc,c;.,-l)=Q^(-.)'^(l + .r'= iC.13) 

Sij, k, c; -1, -1) = f^^^i.fc (C.l^) 

su, k, c; -.-!) = (! + ^Yi-^^f^^r^^^^ ic.i,) 

where, in equation (C.13), the binomial coefEcient (j.) vanishes if j < k. 

Proof. For z = — lorw = — 1, the first two equations follow directly from lemma C.2. Note that and 
^2 have no singularity, in particular at ^ = — 1 or /and w = —1. li w = —1, only the term p = k contributes 

See, e.g., [12] 5.7.1 (7) and 5.10 (3). 
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to equation (C.6), giving equation (C.13). Alternately, setting w = —1 in equation (C.7), the sum over q 
can be performed with lemma C.l. For z = —w — 1, one has zw/ ((1 + z) (1 + w)) = 1, then from the Gauss 
summation theorem 

^ ' ' ^ r{c + j)T{c + k) ^ ' 

ending the proof of equation (C.15). 

Appendix D. Some determinants with gamma functions as entries 

Lemma D. With n a positive integer and c some complex number not zero or a negative integer, let M be 
ann X n matrix which element Mj^^ is equal to T{c + j + k), except possibly for the last row, 

V {Mn-l,k)k=0,...,n-l J 

Then, one has 

^ = r(c + n - 1) y (-1)"-!-^ f " - / , M„_i , {D.2) 

n-1 

det[r(c + j + fc)]^.,^„__^ = n-?!r(c + j). (^.3) 

j=o 

Proof. A determinant is not changed if one adds to any column a linear combination of other columns. 
Replacing for fc = 0, . . . , n — 1 the column Ck of det M by, 

yields a new matrix M' with an equal determinant and which is triangular. Indeed, from definition of M 
and lemma C.l, one has for j = 0, . . . , n — 2 and A; = 0, . . . , n — 1, 

which vanishes for j < k. Then, the determinant is the product of the diagonal elements j\T{c + j) for 
j = 0, . . . ,n — 2 and M^_i Prom equation (D.4), this last element reads 

n-1 / _ 1 \ 1 

M;_i,„_i=r(c + n-l)^(-ir-i-^(^''^ j___M„_i,, {D.6) 

completing the proof of equation (D.2). If moreover, M„_i = T{c + n — 1 + k), using again the equation 
(D.5) for j = A; = n — 1 yields the result (D.3) already published in [15, 14] 

Appendix E. Special cases of Hermite standard bloc orthogonal polynomials 

With the weight functions, x e (—00,00), w := e~^^ and W2 ■= e~^^^, one has: 

See, e.g., [3] 15.1.20 or [13] (1.7.6). 
^9 See [15] equation (A.12) or [14] equation (4.5). 
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n = 1,2... 



1,2,.. 
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{E.l) 
{E.2) 

{E.S) 

(EA) 

{E.5) 



Pi;4 = a;*- Jx^ + i {E.7) 



P0;l 


= X 






Po;2 


= x^- 


1 

4 




Po;3 


= x'- 


3 




Po;4 


= x^- 


3 2 

r + 


3 
16 


Po;5 


= x'- 


5 3 

r + 


15 
16 



i = 2 n = 2,3, .. 
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P2;3 =X^-l^ (E-Q) 

P2;5 = X^-^X^ + lx. (E.IO) 
4 o 

If, instead of w and W2, one considers the following weight functions, 

w^^^x) := w{Vlx) = e-"''' w^^^ := W2{V^x) = e'^""' c> (£.11) 

then, changing x into ^/cx in the integral defining the scalar product ( , )2, the corresponding monic Hermite 
SBO polynomials P^!2 and the orthogonality constant il}.^ read 

pI::!{x) = c-'^-h.n{^^x) = c-(H")^^„ . (E.12) 

In the case i = 1, 'Hermite polynomials under a constraint of zero average' have already been considered 
in [6] and [8], with different choices of weight functions and normalization: 

- in [6] (1) and appendix A, the polynomials An := 2" a;" + 0(a;"~^), n = 1,2, . . ., with w := e~2^ and 
W2 '■= . From equation (E.12) below with c — ^, one gets j4„(x) = 23" Pi;„(a;/\/2), n = 1, 2, . . . ; 

-in [8] (43), the polynomials Q2n := 2""! a;" + 0(a;"-i), n= 1, 2, . . ., with u; := e'^''^ &ndw2 := e"^'. Then, 
one has Q2n = 2~("+^^ A2n, n = 1, 2, Note that with this choice of 102, one has Q2;n = Hn, n = 0, 1, . . .; 

- in [6] section 4, the polynomials Gn ■= CnP\-n, n = 1,2,..., with {Gm , Gn)2 = ^7r/2 5m,n, w := 

and W2 := e^^x^ p^.^^^ equations (2.18) and (2.19), one finds C2„+i = ± 22"-+7^(2n + 1)!, n = 0, 1, ... and 
C2n = ± 22«+i 72^/ V(2n + 1)!, n = 1, 2, . . .; 

- in [8] (44), (45a)-(45d), the polynomials r„ := |c„| Pi-m n = 1, 2, . . ., with c„ as for G„ above, (r^ , r„)2 = 
^/TTi^hm^n-, w :— and W2 ■= e~^^ , and in addition, with a positive coefficient of in r„. 

The 'Hermite polynomials constrained by a zero momentum', considered in [6] section 3 (17), with 
w := e~i^ and 11)2 ■= , is not the case i = 2. Indeed, the one-dimensional subspace £1 considered there, 
and spanned by {x^}, is not 'P2 spanned by {x^,x^}. Nevertheless, as discussed in the conclusion 1-4 (see 
footnote 44 in paper I), from parity arguments, it follows that in this non-standard case, the block orthogonal 
polynomials P„, n = 0, 2, 3, . . . are given by, P2„ = Q2;2n, n = 0, 1, . . . and P2„+i = -P2;2n+i, n = 0,1,.... 
Prom equation (E.12) with c = |, one gets with S„ := 2" a;" -|- 0(a;"~^), n = 0,2,3, .. ., = H2n, n = 
0,1,... and B2n+i(a;) = 23("+5) P2;2„+i(x/\/2), n = 1, 2, . . .. 
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Appendix F. Special cases of Laguerre standard bloc orthogonal polynomials 

With the weight functions, x G [0,cxd), w := e~^x" and W2 ■= e~'^^x", a > —1, one has: 
1 = n = l,2, ... 

Po;i=x-^{a+l) {F.l) 

Po;2=x'^-{a + 2)x+^{a + l){a + 2) (F.2) 

Po;3 =x^- -(a + 3)a;2 + -(a + 2)(a + 3)x- - (q + l)(a + 2)(a + 3) (F.3) 

2 4 o 

^ 3 1 

Po;i = - 2 (a + 4) a;^ + - (a + 3)(a + 4) - - (a + 2)(a + 3)(a + 4) a; 

+ ^(a + l)(a + 2)(a + 3)(a + 4); (F.4) 



n = l,2,... 
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Pi;i =x-{a + l) (F.5) 

A;2 = a;2-i(3a + 5)a;+^(a + l)2 (F.6) 

= x^- {2a + 5) a;2 + i (Sa^ + 19a + 20) a; - ^ (a + l){a^ + 3a + 4) {F.7) 

Pi;4 = - i (5a + 17) a;^ + ^ (Sa^ + 17a + 26) x^ 

- I (7a^ + 48a2 + 125a + 108) x+l{a + lf{a^ + 5a + 12) (F.8) 

8 8 

P1.5 = x^ - {3a + 13) + ^ {To? + 53a + 106) -\ {^o? + 39a2 + 137a + 162) 
+ ^ (9a'' + 98a3 + 447a2 + 886a + 648) x 

- ^ (a + l)(a^ + lOa^ + 47a2 + 86a + 72) ; (F.9) 
i = 2 n = 2,3, ... 

P2;2 = - 2(a + 2) a; + (a + l)(a + 2) (F.IO) 

-P2;3 = - ^(5a + 13) + 2(a + 2)^ x - ^(a + l)2(a + 2) (P.U) 

P2.4 = x^ - (3a + 10) x^ + ^(13a2 + 71a + 102) x^ 

- ^(a + 2)(3a^ + 13a + 18) x + ^(a + l)(a + 2)(a^ + 3a + 6) {F.\2) 
P2;5 = x^ - ^(7a + 29) x^ + ^(19a2 + 135a + 254) x^ 

- l{25a^ + 222a^ + 719a + 810) x^ + \{a + 2 f{2a^ + 11a + 27) x 
8 2 

In the case i = 1, 'modification of Laguerre polynomials by a constraint of zero average' have already 
been considered in [7] section 2, with different choices of weight functions and normalization: the monic 
polynomials G'|^, n = 1,2,..., with a := d — 1, w := e~^^ x'^~^ and W2 := x'^^^ . Setting c = ^ in 

equation (F.16) below yields: G'^{x) = 2" Pi;„(x/2)|a=d_i, n = 1,2, Note that with this choice of W2, 

Q2;n = L^n \ n = 0, 1, . . .. 
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-l{a + lf{a + 2){a'^ + 5a + 18) {F.13) 
o 

P2;6 = - 4(a + 5)x^ + ^(13a^ + 115a + 268) - ^(Ha^ + 127^2 + 526a + 752) 

+ ]^(41a!^ + 538a3 + 2947a^ + 7418a + 7056) 

-l{a + 2) (5a* + 58a^ + 319a^ + 770a + 720) x 
8 

+ -^(a + l)(a + 2)(a'* + lOa^ + 59a2 + 122a + 144) . (F.14) 
If, instead of w and W2, one considers the following weight functions, 

^(''^(a;) := c-"m;(cx) = e'^'^x" w^^'^^x) := c-°'w2{cx) = e'^^^a;" c> {F.15) 



then, changing x into cx in the integral defining the scalar product ( , )2, the corresponding monic Laguerre 
SBO polynomials Pj^.^^ and the orthogonality constant H'j^^^ read 

^^(a.) = c-^hAcx) flg = c-(«+i+^")H,;„ . (i^l6) 
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